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Abstract— We review the use of graphene to develop 
reconfigurable, miniaturized, and efficient terahertz 
(THz) antennas and associated feeding networks, and 
attempt to identify current research trends and mid- and 
long-term challenges and prospects. We first discuss the 
state of the art in resonant, leaky-wave and reflectarray 
antennas, providing a critical assessment of their 
performance, limitations, and main challenges that 
remain to be addressed. Next, we examine different 
integrated feeding networks, including components such 
as switches, filters, and phase shifters, and we clarify the 
impact that graphene's intrinsic spatial dispersion may 
have in their performance. Our outlook clearly describes 
how graphene can bring exotic functionalities to all these 
devices, including quasi real-time reconfiguration 
capabilities and magnet-less non-reciprocal responses. 
Some exciting applications of THz antennas are then 
presented and discussed, including transceivers, 
biosensors, and first experimental realizations of detectors 
and modulators. We conclude by outlining our vision for 
the promising future of graphene-based THz antennas. 


Index Terms—Terahertz, graphene, dipoles, leaky-wave 
antennas, reflectarrays, reconfigurability, non- 
reciprocity, sensors. 


I. INTRODUCTION 


Advanced electromagnetic antennas operating at micro- and 
millimeter wave frequencies have become indispensable in 
today’s information society [1], [2]. They are massively 
employed in a wide range of applications including 
multimedia broadcasting, mobile or satellite communications, 
radars, environmental sensing, or medical systems [3]-[14]. 
This plethora of applications has naturally triggered enormous 
research and technological progress in the field in the past 
decades. However, the exponential growth of the information 
society is continuously imposing very challenging — and 
sometimes conflicting — technological demands, such as 
higher communication data rates, higher number of 
functionalities in a single, miniaturized device, future 5G 
heterogeneous networks, fully integrated and miniaturized 
sensing solutions, and the efficient use of the spectrum and 
energy resources, to name a few. 


One straightforward solution to fulfill such demands would 
be to further increase the operation frequency of current 
antennas, taking advantage the terahertz (THz) spectrum [15]— 


[18], i.e., electromagnetic (EM) waves with frequencies 
ranging from 0.1 to 10 THz [16], [19]. In fact, this band 
presents unique opportunities for advanced applications [16], 
[19], including wireless communications with hundreds of 
gigabit per second or even terabit-per-second data rates, high- 
speed miniaturized processing systems, sensing the presence 
of specific chemical or biological processes, ultra-high 
resolution imaging systems for inspection, screening, and 
non-invasive identification of materials in real-time by 
spectroscopy, study of cosmic background radiation, or 
detecting gases in the Earth and in the atmosphere of other 
planets, among many others [16]-[22]. Despite recent 
promising advances [21], these applications are not currently 
being fully-exploited by the information society due to the 
immature state of THz technology in terms of antennas, 
detectors, sources, and basic components. In fact, this 
frequency region is known as the “THz gap”, as it occupies a 
technology gap between the well-developed areas of 
electronics and photonics. 


In a related context, graphene, a two-dimensional carbon 
material, has recently triggered very intense and 
multidisciplinary research efforts due to its outstanding 
electromagnetic, mechanical, electrical, and thermal 
properties [23]-[32]. Graphene supports the propagation of 
surface plasmon polaritons (SPPs) able to exhibit strong wave 
localization, moderate loss, and the exceptional property of 
being tunable through electrical/magnetic bias or chemical 
doping. Importantly, since such plasmonic response appears 
in the terahertz and infrared frequency bands, graphene has 
quickly positioned itself as a very promising platform for THz 
transceivers and optoelectronic systems. Even though SPPs 
can be supported by composite structures at THz [33], 
graphene is yet the only pristine material able to provide such 
plasmonic response, thus opening exciting and unexpected 
prospects for wave manipulation and radiation in this 
frequency band. Graphene’s unique properties have been 
exploited to propose miniaturized and reconfigurable resonant 
[34]-[40], leaky-wave [41]-[46], and reflectarray antennas 
[47]-{[50] with unrivaled radiation efficiency and 
functionalities at THz. These antennas resonate at frequencies 
much lower than their metallic counterparts and they fully 
exploit the reconfigurable nature of graphene’s conductivity 
[51]. In addition, all fundamental building blocks of front-end 
transceiver architectures can be realized using graphene 
plasmonics, including switches, phase-shifters, low/band pass 
filters, or modulators. Non-reciprocal propagation and giant 
faraday rotation have also been demonstrated when graphene 
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is biased by magnetostatic fields [52]—[54], but the bulky 
magnets required with this approach are antithetical to the 
desirable miniaturization of graphene plasmonics. Even more 
excitingly, graphene’s unique properties have led to the 
development of magnet-free non-reciprocal guided devices 
and antennas based on the spatiotemporal modulation of its 
conductivity, allowing to break the fundamental principle of 
time-reversal symmetry (reciprocity) without ferromagnetic 
materials or magneto-optical effects [43], [46], [55]—[62]. The 
adequate integration of all these components and 
functionalities into a single, miniaturized, device will 
significantly push the boundaries of current terahertz 
technology, paving the way towards silicon-compatible THz 
communication and sensing systems with unprecedented 
performance and responses. 

It is important to mention that, to date, most graphene-based 
antennas, devices, and detectors have only been proposed, 
designed, and analyzed theoretically, and their 
electromagnetic response have been predicted using 
numerical simulations. However, we do highlight that recent 
experiments have fully verified the propagation of confined 
and tunable SPPs in graphene [63]-[66], which is the physical 
phenomenon that all antennas and components described here 
rely on. Earlier experiments yielded graphene of very poor 
quality that was not viable for any electrodynamic application 
other than to serve as a tunable resistor, severely limiting the 
practical feasibility of most research. The past few years, 
however, have seen immense progress in this regard, thanks 
to extensive efforts by dozens of groups worldwide to improve 
fabrication techniques. For instance, mobility improvements 
of several orders of magnitude have been reported for 
graphene structures sandwiched between hexagonal boron 
nitride layer [63], [67]-[69], which are well above the quality 
requirements of all antennas and devices discussed in this 
review. More recently, and despite the challenging fabrication 
and measurement techniques involved, the response of some 
exciting THz detectors [70], [71] and modulators [72]-[78] 
have already been experimentally confirmed, clearly 
outperforming previous state of the art and showing excellent 
agreement with numerical expectations. These pioneering 
experimental works have just barely unveiled part of the 
immense potential of THz plasmonics using graphene or even 
other two-dimensional materials such as black phosphorus 
[79]-[83] or 2D chalcogenides and oxides [84]-[89], 
confirming that this exciting field is still in its infancy and that 
the future of graphene and 2D materials in antenna research 
looks brighter than many ever dared to imagine during its 
early steps 

In the following, we summarize and contextualize the 
evolution of THz graphene-based antenna research since its 
beginnings merely 6 years ago, and attempt to identify current 
research trends and mid- and long-term challenges and 
prospects. We review the state of the art in resonant antennas, 
leaky wave antennas and reflectarrays, and conclude with an 
overview of their feeding networks and the most exciting 
applications envisioned for graphene antennas and THz 
systems. 
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Fig. 1. Complex conductivity of graphene evaluated using the Kubo 
formula for several chemical potentials p, and ratio between its real and 
imaginary parts. Graphene relaxation time is T = 0.5 ps. 


IL. ABOUT GRAPHENE CONDUCTIVITY 


Graphene’s conductivity directly depends on its unique 
band structure [26], [90]-[92] and a set of parameters such as 
chemical doping, Fermi energy (chemical potential), electron 
mobility, relaxation time, operation frequency or the 
environmental temperature. We stress that these parameters 
are not independent, but might be related to each other [26], 
[30], [90], [93], [94]. Using a semiclassical, local model and 
in the absence of magnetic bias, graphene's conductivity can 
be modeled at THz frequencies just using the intraband term 
of the Kubo formula [25], [95] 

"ES jin {2 f + cosh ( a Jl (1) 
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where e is the electron charge, T is the electron relaxation time 
in graphene, kg is Boltzmann's constant, T is temperature, h 
is the reduced Planck's constant, w is the angular frequency, 
and p, is graphene's chemical potential. This expression is 
valid if hw < 2|u,|, which holds at THz frequencies for 
moderate uç. At higher frequencies, interband contributions of 
graphene conductivity becomes significant and must be 
considered [26]. An example of graphene’s dispersive 
conductivity is shown in Fig. 1. An exhaustive overview of 
such conductivity and associated phenomena is out of the 
scope of this paper. Instead, we refer the reader to some of the 
excellent literature available on this topic [25], [91], [92], 
[95], [96]. For the sake of completeness, we will briefly 
elaborate on some of its exciting properties that have direct 
application in antenna research. 


From an electromagnetic perspective, graphene is 
frequency-dependent and behaves as a surface resistor at 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


micro and millimeter wave frequencies. However, its 
conductivity acquires a different response at terahertz and 
infrared frequencies because its intrinsic kinetic inductance 
(associated to the imaginary part of the conductivity) plays the 
role of negative real permittivity in a bulk material. 
Consequently, and similarly to noble metals at optics, 
graphene supports the propagation of SPPs at these 
frequencies, which has recently triggered the development of 
antennas and guided devices. Specifically, in a free space 
standing graphene sheet of conductivity o, the dispersion 
relation and characteristic impedance of the supported 
transverse magnetic (TM) plasmons can be computed as [25], 
[96], [97] 


2\2 
spp — ks ]1 — (2) Q) 
_ _Kspp 
"e sur (3) 


where ky and 7 are the free-space wavenumber and 
impedance, € is the permittivity of vacuum, and égr¢ is the 
effective permittivity of the surrounding media. Even though 
this simple model assumes laterally infinite sheets, it is 
nonetheless an extremely useful estimate of the dominant 
plasmon wavelength and loss in realistic structures like 
nanoribbons and patches. 


Undoubtedly, one of the most exciting properties of 
graphene is the possibility to manipulate its complex 
conductivity. This can easily be done exploiting graphene’s 
field effect by simply applying an external electrostatic field 
bias perpendicular to the sample [98]. The bias injects 
electrons or holes on the sheet, thus modifying graphene’s 
chemical potential (see Fig. 1). Recent experiments have 
confirmed that this tunability can be very significant and fast 
— up hundreds of gigahertz [73], [99] -, thus enabling 
reconfiguration capabilities in all graphene-based antennas 
and components. In addition, when a magnetostatic bias is 
applied to graphene, its band structure is not continuous 
anymore but discretized in the so-called Landau levels [100]- 
[109]. Graphene’s conductivity becomes then a gyrotropic 
tensor and leads to exciting non-reciprocal responses such as 
Faraday rotation or isolators [110], [111]. Optical pumping 
has also been used to modify graphene’s properties, enabling 
THz lasing and operation as reconfigurable gain media [112]- 
[115]. It is also worth mentioning the non-local response of 
graphene [96], [116]-[118], which becomes important when 
the supported ultra-confined SPPs presents variations at 
speeds comparable to the carriers on the materials. This 
phenomenon must be considered when designing guided 
devices, as described later, since it can significantly alter their 
expected response. 


III. RESONANT ANTENNAS 


Resonant antennas are ubiquitous in daily life, so it is only 
natural that they were the first kind of graphene antennas 
theoretically investigated. Pioneering works by Jornet and co- 
workers demonstrated that subwavelength graphene patches 
can resonate at THz frequencies, making them very attractive 
antennas for this elusive frequency band [34], [35], [119]. In 
their initial studies about graphene nano-ribbon (GNR) 
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Fig. 2. Graphene-based patch as a THz scatterer [34]. (a) 3D 
schematic of the structure [39]. (b) Resonant frequency of the patch 
versus its geometrical features [39]. Graphene chemical potential and 
relaxation time are 0.0 eV and 0.1 ps, respectively. Reprinted from [39] 
with permission from Elsevier. 
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Fig. 3. Single layer graphene-based planar dipole [37]. (a) 3D schematic 
of the structure. A THz photomixer — red region with width G — is 
employed to excite the antenna. Inset shows the E-plane radiation pattern. 
(b) Input impedance of the graphene dipoles described in Table 1. © 
Reprinted from [40], with the permission of AIP publishing. 


dipoles, they did not yet consider the direct 
detection/excitation of | SPPs for operation in 
reception/transmission, studying instead the extinction cross 
section of the structure to determine its resonant frequency. 
The configuration they considered is depicted in Fig. 2, and is 
comprised of a graphene ribbon transferred onto a substrate of 
thickness D. Fig. 2b shows the resonant frequency versus the 
ribbon width W and length L. The ability of graphene patches 
and strips to resonate in the THz band arises from the 
unusually high kinetic inductance of graphene, which enables 
short-wavelength SPPs in this frequency range [120]-[122]. 
The associated resonant frequencies can easily be 
approximated by considering the propagation of SPPs along 
patches of finite length and then computing their Fabry-Perot 
resonances [39]. Examples of similar analysis using common 
transmission line techniques will be further described below. 
Once these resonant frequencies are available, the effective 
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Fig. 4. Double layer graphene-based planar dipole [38] . (a) 3D 
schematic (top) and lateral view (bottom) of the structure. (b) E-plane 
radiation pattern with (solid blue line) and without (dashed red line) 
including a silicon lens to the antenna. The lens has a height of 572 um, 
a radius of 547 um and is made of Silicon modelled with £,. — 11.66 
and tanó — 0.001. The dielectric thickness is 72160 jm. (c) Input 
impedance of the graphene antenna upon variation of the chemical 
potential. Radiation (d) and total (e) efficiency — considering a 10 kQ 
photomixer as a source — at resonance versus chemical potential. The 
antenna has a total length Z=11 wm and a width W=7 pm and is 
composed of two stacked graphene patches [76] separated by a h=100 
nm thick AbO; layer (£, — 9 and tanó — 0.01). © Reprinted from 
[38], with the permission of AIP publishing. 








Hc(eV) | L(um) | W(um) | W.P. | f(THz) | L/Ag 
0.13 17 10 L 1.023 | 0.06 
H 1.35 0.08 

0.25 25 20 L 1.472 | 0.09 
H 1.534 | 0.12 


























Table 1. Characteristics of the single layer graphene-based planar 
dipoles studied in Fig. 3. 


polarizability of the subwavelength patch can be retrieved, 
thus allowing to easily obtain the absorption and scattering 
cross sections of the structure. 
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Fig. 5. Circuit model of graphene-based plasmonic dipoles [124]. © 
2014 IEFE. Reprinted, with permission, from [124]. 


Soon after this study, Tamagnone and colleagues proposed 
the first true graphene antenna, in the sense that it bridges the 
gap between SPPs propagating on the sheets and far field 
radiation [37]. The geometry considered, depicted in Fig. 3, is 
similar to Fig. 2 and is based on the same operation principle. 
The key contribution is the addition of a port able to efficiently 
excite SPPs over graphene (or detect them in reception, from 
reciprocity). This allowed to perform the first quantitative 
study of input impedance and radiation efficiency of GNR 
dipole antennas, demonstrating unprecedented tunability via 
electrostatic bias and high values of input impedance suitable 
for THz sources such as photomixers. As an example, the 
input impedance of a graphene dipole antenna tuned to He = 
0.13 and 0.25 eV is shown in Fig. 3b. The frequency points 
suitable for operation as an antenna (i.e., real input 
impedance) are marked in the plot as WPL and WPH (working 
points low and high, respectively), corresponding to the 
conditions L ~ 0.5Agpp, and L % Àçpp where Àspp is the 
wavelength of the SPP mode supported by the graphene strip. 


As shown in Table 1, the dimensions of this antenna (=) 
0 


remain deeply subwavelength. These promising findings 
encouraged the systematic study of graphene at THz from the 
unexplored perspective of antenna engineering. 


This graphene patch antenna was refined shortly after by 
considering a double layer graphene stack and including a 
dielectric lens to increase its directivity [38] (see Fig. 4). This 
way, the advantages of wide tunability and real input 
impedance in subwavelength plasmonic devices could also be 
exploited in applications requiring directive antennas. 
Moreover, they demonstrate excellent impedance stability 
under wide reconfiguration of the resonant frequency, as 
shown in Fig. 4c. A more comprehensive study on the 
radiation and total efficiency of the proposed antenna as a 
function of chemical potential is summarized in Fig. 4d-e. We 
do highlight that the computed radiation efficiency is indeed 
high, clearly outperforming other antenna of similar size 
operating at THz [123]. Besides, it is interesting to note that 
the efficiency, unlike the input impedance, is not stable upon 
reconfiguration. This is easily explained in terms of electrical 
size of the antenna. Increasing the chemical potential pu, 
increases the plasmon wavelength, upshifting the resonant 
frequency and increasing the electrical size of the antenna in 
terms of free-space wavelength. This is a common theme in 
all graphene antennas: extremely miniaturized devices can be 
designed by exploiting the high inductance of graphene for 
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small u,, but it usually comes at the cost of radiation 
efficiency. A rarely mentioned downside of using large pu, is 
that the properties of graphene SPPs become less sensitive to 
bias voltages, negatively affecting electronic 
reconfigurability. 


The design and analysis of graphene-based dipoles 
discussed so far relied heavily on full-wave simulations to 
compute the input impedance and the resonant frequency. In 
2014, a simple and elegant circuit model for this type of 
antenna was proposed [124]. This model, depicted in Fig. 5, 
introduces a capacitance between the two dipole arms and a 
reactance at their termination to capture the effect of fringing 
fields. The characteristic of the transmission line — which 
represents SPPs propagation through the graphene patch — can 
be calculated exactly from numerical eigensolvers, or 
approximately by Eqs. (1)-(3) for the case of laterally infinite 
sheets. With these elements, the antenna resonant frequency, 
input impedance, and radiation/total efficiency can be 
calculated using straightforward transmission line theory. 


To further increase the radiation efficiency of these 
antennas, the concept of graphene coated wires (GCW, see 
Fig. 6) antennas was recently introduced [40]. Authors 
analytically developed the dispersion relation of multiple 
variations of this cylindrical waveguide configuration 
supporting electrical and magnetic biasing, and showed that 
the absence of edges substantially reduces the propagation 
loss of SPPs compared to strip-based configurations [40], 
[125]-[128]. Moreover, this class of waveguides can be 
designed to be immune to the dielectric loss of the wire 
material, since all power travels outside [125]. These findings 
have important implications in the design of optimal THz 
interconnects and reconfigurable waveguide devices, but also 
in the design of resonant antennas. While the patch-based 
antennas have the inherent advantage of planarity, GCW 
dipole antennas were shown to present significantly higher 
radiation efficiency than similar planar implementations. 
Table 2 presents a performance summary for two comparable 
planar and GCW resonant antennas. GCWs presents 
significantly better performance, attributed to the their higher 
SPPs quality factor and increased robustness to dielectric loss 
and possible trapping of radiation by the substrate. Fig. 6c-d 
depict the input impedance of similar GCW- and patch-based 
dipoles, demonstrating an extended reconfiguration range for 
the GCW. Fig. 6e-f show the input impedance and radiation 
resistance [1] computed for a wide range of frequencies and 
Lc. The maxima and minima of radiation resistance in Fig. 6f 
roughly correspond to the odd and even multiples of the 
plasmon wavelength along the dipole, respectively. Like all 
resonant dipoles based on deeply confined plasmons, the 
radiation nulls appear because the contributions of adjacent 
subwavelength current units tend to cancel out in the far-field, 
which is not the case in conventional dipole antennas [40], 
[129], [130]. We note that a circuit model similar to Fig. 5 was 
used to design and characterize this class of resonant antennas 
(thin traces in Fig. 6c), further demonstrating the general 
validity of traditional microwave concepts in modeling 
graphene antennas. 


We have presented an overview of the most fundamental 
graphene-based resonant antennas, but many groups 
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Fig. 6. Reconfigurable plasmonic antennas based on graphene tubes 
[40]. (a) 3D schematic of a graphene-based cylindrical waveguide. Insets 
illustrate the lateral view of a single (left) and dual (right) layer graphene 
cylindrical waveguides. (b) Operation of the device as an antenna. (c) 
Input impedance of a double-layer graphene cylindrical dipole upon 
reconfiguration. Thick traces computed with CST MWS, and thin traces 
computed with an equivalent circuit (see Fig. 5 and [40]). Parameters are 
L=26 pm, a=450 nm, b=500, €,1 = £,4 — 3.8. (d) Input impedance of 
a similar-sized double-plane graphene planar antenna (see Fig. 3 and 
[38]). Parameters are W=1 um, permittivity and thickness of substrate 
€, = 3.8 and t=1 pm, and a distance between graphene layers h=50 nm. 
Real part of the input impedance (e) and radiation resistance (b) of the 
plasmonic dipole defined in panel (c) versus frequency and chemical 
potential. © 2015 IEEE. Reprinted with permission, from [40]. 


























GCW Strip 
Hc(eV) | Vo(V) Trad | "hot | Vo(V) | Mraa | Not 
0.05 0.46 12 1.1 0.43 0.16 | 0.12 
0.1 1.84 4 3.7 1.75 0.5 0.4 
0.2 7.3 14 12.5 7.0 2 1.4 
0.3 16.5 26 24 15.7 3.8 3.5 





























Table 2. Bias voltages (Vo in Fig. 6a) and efficiencies (96) of the 
graphene cylinder- and strip-based antennas described in Fig. 6 and Fig. 


worldwide have pursued ways to improve their performance 
or to achieve exotic responses. The variations are too 
numerous to describe in detail here [131]-[138], but include 
such interesting concepts as enhancing the radiation of 
quantum systems [139], an efficient integration of a 
photomixer in the graphene antenna [140], or non-reciprocal 
behavior based of magnetic bias with well-defined 
fundamental bounds [141], [142]. 
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Fig. 7. Metal-graphene hybrid antenna [143]. (a) Top (left) and lateral 
(right) view of the structure. (b) Total radiation efficiency — considering 
a 10 kQ photomixer as a source — of the hybrid antenna upon 
reconfiguration. © 2013 IEEE. Reprinted with permission, from [143]. 
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Fig. 8. Reconfigurable MIMO system based on graphene patch 
antennas [147]. (a) Schematic of the array. Each element is composed 
of a graphene layer with width W=2 um and length L=5 um. (b)-(d) 
Radiation pattern of the antenna at 1 THz. (b) Elements 1, 3, and 4 
are in a low-resistance state. (c) Elements 1, 3, and 5 are in a low- 
resistance state. (d) Elements 2, 3, and 5 are in low-resistance state. 
Graphene chemical potential and relaxation time are uw, = 1.0 eV and 
T = 0.5 ps, respectively. © 2014 IEEE. Reprinted with permission, 
from [147]. 
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A different approach to exploit the tunability of graphene’s 
conductivity in the design of THz antennas consists on using 
it as a perturbation in other structures with potentially higher 
radiation efficiency, typically involving metals [143]. The 
goal is to use graphene to electrically tune the resonant 
frequency, but not as the (main) radiating element. This 
concept is illustrated in Fig. 7 [137]. The graphene area (red 
color) located at the antenna center has the same geometry as 
the previously discussed planar dipoles, and is connected to 
the source. The yellow area represents metal, designed to 
maximize radiation efficiency while efficiently coupling to 
the graphene strips. We highlight that graphene’s main role 
here is to serve as a matching network between the photomixer 
source — which presents a very high input impedance — and the 
metallic structure — which presents a high radiation efficiency. 
By tuning graphene’s conductivity like in previous planar 
dipoles, the input impedance of the whole structure can be 
controlled in a wide range, tuning the total efficiency of the 
antenna — the radiation efficiency is barely affected since 
graphene is not the main radiator. Fig. 7b shows the antenna 


total efficiency, demonstrating the wide tunability and 
increased in efficiency compared to the graphene-only dipole 
of Fig. 3-4, obtained at the cost of significantly increased size. 
This example once again illustrates the recurring trade-off 
between miniaturization and efficiency in graphene antennas. 
We stress that many other works have also considered the use 
of graphene as a parasitic element, for instance using it as a 
tunable substrate [144], [145] or as a high-impedance surface 
[146]. 

All graphene resonant antennas discussed in this section can 
be used as part of more complex systems of arrays. As an 
example, Fig. 8 illustrates its use in a multiple input, multiple 
output (MIMO) configuration [147]. The system is comprised 
of several graphene nano-patch antennas whose resonant 
frequencies can be independently controlled. This 
configuration permits to easily achieve efficient and tunable 
MIMO responses at THz while keeping a very compact 
geometry and integrated, subwavelength elements. Fig. 8b-d 
illustrate the capabilities of the system to produce different 
radiation patters by changing the conductivity of individual 
patches in an array of just a few elements. 


IV. LEAKY-WAVE ANTENNAS 


Leaky-wave antennas (LWAs) have been a very active topic 
of antenna research in the past decades due to their interesting 
beam-scanning properties [148]-[159]. They are based on the 
energy leaking towards free-space from an electromagnetic 
wave that is propagating through a waveguide and whose 
phase velocity is greater than the speed of light. Even though 
most LWAs have been designed to operate at micro and 
millimeter waves, their underlying physics is fundamental to 
explain exotic physical concepts such as extraordinary 
transmission [160], [161], electromagnetically induced 
transparency [162], or Cherenkov radiation [163]. Graphene 
is particularly suited for manipulating leaky waves at THz and 
infrared frequencies, allowing to implement reconfigurable 
structures and, more importantly, to achieve unprecedented 
magnet-less non-reciprocal responses. 


One common and powerful type of leaky-wave antennas is 
based on the sinusoidal modulation of reactance surfaces. This 
concept was proposed by Prof. Oliner in the 50s [164] and has 
recently received significant attention at microwaves [165]- 
[168]. Essentially, the periodic modulation allows to express 
the fields along and on top the surface as an infinite sum of 
Floquet harmonics, some of which may lie within the light 
cone (k < ko) and become leaky waves, causing the structure 
to operate as an antenna. The main parameters that control this 
type of antennas are the average reactance X,, the amplitude 
of the modulation M, and the modulation period p. The 
reactance should be inductive to support confined TM waves 
that strongly interact with the modulated surface — its specific 
value determines the wave propagation constant and 
confinement. The modulation amplitude and period control 
the radiation rate and angle, respectively. The effective 
modulated surface reactance (assuming modulation and 
propagation along y) can thus be expressed as 


X'; = X,[1 + Msin(2ry/p)]. (4) 
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Fig. 9. Graphene leaky-wave antennas based on a sinusoidally 
modulated impedance [41], [42]. (a) 3D schematic of the first 
implementation, which is composed of a graphene sheet transferred onto 
a back-metallized substrate and several polysilicon DC gating pads 
located beneath it [42]. Spatial modulation is achieved through electrical 
gating. (b) Radiation diagram for different lengths of the electrically- 
induced structure periodic modulation [A, = N- (s + g)]. Operation 
frequency is 2 THz, g=0.2 jm, and the applied voltage is in the V pç = 
6.5— 45 V range. Top (bottom) diagrams have been generated 
considering an antenna transferred over a substrate with €, = 3.8 (1.8) 
and gating pads with s=4.8 (8.8) wm. (c) 3D schematic of the second 
implementation, which achieves the spatial modulation through the 
modulation of the graphene strip width [41]. (d) Radiation diagram at 
1.5 THz for different chemical potentials. Results have been computed 
using full-wave simulations (solid lines) and leaky-wave theory (dashed 
lines). Graphene relaxation time is T = 1.0 ps. Panels (a), (b): © 2014 
IEEE, reprinted with permission, from [42]. Panels (c), (d): © 2013 
OSA. Reprinted with permission, [41]. 


Graphene is the perfect platform to implement this concept at 
THz frequencies, since it provides an inductive response, 
given by Z, =1/o =R;+jXs, where R; and X; are the 
surface resistance and reactance of the sheet, respectively, that 
can be manipulated at will by exploiting graphene’s field 
effect [98]. 


The pioneering first graphene-based LWA, proposed by 
Esquius-Morote and co-workers in 2014, is illustrated in Fig. 
9a-b [42]. The structure is composed of a graphene sheet 
transferred onto a back-metallized substrate and a set of 
polysilicon DC gating pads located beneath it. By applying 
independent voltages to each pad, arbitrary modulations of the 
reactance profiles can be implemented over the sheet, thus 
allowing the propagation of leaky-waves. Importantly, the 
pointing angle 6) and leakage rate a,,q of the radiated beam 
can be dynamically controlled by changing the features of the 
modulation. The voltage required in each pad depends on the 
substrate capacitance and the average reactance (average Lt), 
and can easily be found analytically [42]. The wavenumber ky 
of the fundamental harmonic along the modulated graphene 
sheet can be written as 


ky = Bspp + AMBspp — J(&spp t araa) (5) 


where [spp is the SPP propagation constant on the 
unmodulated sheet (M=0), ABspp is a small variation in the 
wavenumber due to the modulation (usually negligible), aspp 
the dissipative attenuation constant, and a, qq the attenuation 
constant due to leaky radiation. The space harmonic typically 
used for radiation in periodic leaky wave antennas is the -1 
harmonic, with propagation constant 6_, = k, — 2m/p. The 
pointing angle of the radiated beam is then 


. ~ acin [BSPP _ 20 
Oo = asin(B_,/ko) ~% asin (Ere — =) (6) 

Fig. 9b demonstrates outstanding capabilities of the antenna 
to switch from backward to forward radiation and quasi- 
broadside by simply changing the gate voltages. We do note 
that maximum radiation pointing exactly at broadside is not 
possible due to the open stopband of the periodic structure 
[148], [156]. The relatively lossy nature of graphene plasmons 
makes it challenging to simultaneously achieve highly 
directive beams and high-radiation efficiency, and one might 
decide to favor over the other through the choice of the 
modulation amplitude: small M for high effective aperture and 
thus higher directivity, or large M for increased efficiency. We 
do stress, however, that this antenna, with realistic parameters 
of tT = 1 ps and M = 0.35, has an impressive (for the THz 
band) radiation efficiency of roughly 11% while 
simultaneously providing beamscanning features. 


A different implementation of graphene-based LWAs at 
THz was also proposed by same authors [41]. The antenna is 
shown in Fig. 9c-d and instead of gate pads uses a periodic 
modulation of the strip width to excite leaky spatial 
harmonics. To this purpose, authors take advantage of the 
large sensitivity of plasmons propagation constant with 
respect to the strip width. For ease of design and calculation, 
this structure can be modelled as an effective modulation of 
the reactance by simply mapping the width-dependence of the 
plasmon dispersion to effective changes on graphene’s 
conductivity o = 1/Z, . Such analysis can efficiently be 
performed using the analytical scaling law of SPPs versus 
strips width introduced in [97], which was applied to realize 
quick designs and analysis of this type of antennas in [41]. 
This implementation is significantly simpler to fabricate than 
the previous one, since it avoids the use of gating pads beneath 
the graphene sheet, while still allowing a great control over 
the radiated beam direction at a fixed frequency by applying a 
single DC bias voltage (see Fig. 9d). There is, however, less 
control over the modulation amplitude M and operation 
frequency than in the previous configuration, since they are 
determined by the width modulation at the fabrication stage. 


Other approaches to implement LWAs include borrowing 
concepts from the microwave domain, such as composite 
right-left handed transmission lines and antennas [152], [153]. 
Specifically, [45] studied the viability of implementing such 
techniques at THz frequencies using graphene plasmonics. 
Even though the preliminary results are somewhat promising, 
the required structures are quite challenging to fabricate and 
need very high-quality graphene samples, which might delay 
its use in practice. 


Elastic vibrations based on flexural (mechanical) waves 
have also been proposed as an alternative method to impart 
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Fig. 10. Graphene leaky-wave antenna based on a sinusoidally 
modulated impedance achieved through elastic vibrations from flexural 
waves [44]. Inset (top right) shows the antenna radiation diagram at 30 
THz versus different flexural frequencies. Graphene chemical potential 
and mobility are 0.5 eV and 20,000 cm?/V's, respectively. © 2014 IOP 
Publishing. Reprinted with permission, from [44]. 


the periodic modulation required to excite leaky waves [44], 
fully exploiting the multidisciplinary nature of graphene. This 
approach, illustrated in Fig. 10, relies on the fact that a flexural 
wave travelling on graphene can be accurately modeled as a 
static grating with period p (A in the figure). By changing the 
frequency of the flexural waves with a biharmonic source, the 
radiated beam can be dynamically controlled. This antenna 
has the main advantage of providing any desired modulation 
reactance profile, thanks to the flexibility of mechanical 
waves. In fact, its analog nature inherently avoids the 
discretization of the reactance profile using gate electrodes or 
its fixing during fabrication. Consequently, it can implement 
almost any desired radiation pattern and provide 
beamscanning functionalities. However, the structure is not 
fully planar and might be difficult to fabricate, since it requires 
free-standing graphene able to vibrate. Moreover, the 
adequate combination of miniaturized plasmonic devices and 
mechanical components may be challenging, especially in 
practical scenarios. 


Another possibility to implement graphene-based LWAs is 
to use graphene as a perturbation element, as previously 
discussed for the case of resonant antennas. The motivation is, 
once again, to exploit the advantages of different technologies 
to maximize performance or reduce the restrictive dependence 
on high-quality graphene. In these designs, graphene is almost 
invariably used to provide quasi real-time electronic 
reconfigurability that would be almost impossible to achieve 
with conventional materials. This is the case of LWAs based 
on high-impedance surfaces [146], [169], which achieved 
enhanced beamforming and tunable features thanks to an 
adequate use of graphene. A recent example of an alternative 
approach is illustrated in Fig. 11, where graphene is included 
in a substrate-superstrate leaky wave antenna to allow fixed- 
frequency scanning [170], [171]. Here, the position of the 
graphene sheet is optimized such that it significantly affects 
the propagation constant of the supported leaky modes. After 
a careful optimization, the structure exhibits strong 
reconfigurability, minimizes the leakage rate ay to maximize 
directivity, and generates tunable radiation with nearly 
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Fig. 11. Substrate-superstrate graphene-based leaky-wave antenna 
[170], [171]. (a) Cross section of the structure. A graphene sheet is 
located in the substrate, and a horizontal magnetic dipole is employed to 
excite the antenna from the ground plane. (b) Radiation patterns (E- 
plane) at 1.14 THz upon reconfiguration. Parameters: €,, = 3.8 (SiO2), 
£,2 — 25 (HfO;), xo/h, — 0.81, and graphene relaxation time is 0.5 
ps. © 2016 IEEE. Reprinted with permission, from [171]. 


constant beamwidth. Fig. 11b demonstrates the 
reconfiguration capabilities of this antenna, with up to 45 deg 
of beam scanning by tuning the chemical potential from 1 eV 
to 0.3 eV using a single electrode. 





In a different context, a recent and very promising trend in 
electromagnetics and in the antenna community is the pursuit 
of non-reciprocal responses without ferromagnetic materials 
or magneto-optic effects. This functionality traditionally 
requires bulky and expensive magnets that severely limit the 
miniaturization of devices and their integrability with other 
technologies. This is especially important in the case of THz 
graphene plasmonics, where the resulting components are 
significantly smaller than the silicon-incompatible required 
magnets. A new paradigm has recently been proposed to break 
time-reversal symmetry by modulating the properties of a 
structure (permittivity) in space and time. This approach has 
been successfully applied across the EM spectrum [56], [58], 
[59], [61], acoustics [55], [60], and even mechanics [172], 
allowing, for instance, the development of exciting circulators 
with state of the art performance that are compatible with 
integrated circuits [58]. A transition towards a magnet-less 
non-reciprocal paradigm will enable cheap and compact 
versions of essential devices such as isolators, circulators, and 
non-reciprocal antennas. 


In the exciting work of Correas-Serrano [43] such concepts 
are translated into graphene plasmonics at THz by adequately 
exploiting the excellent reconfigurability of the material 
through its field effect. Specifically, this approach enables the 
development of THz antennas with different transmission and 
reception properties, as illustrated in Fig. 12. Even more 
interestingly, under time reversal, the SPPs 
transmitted/received by the device oscillate at different 
frequencies. The combination of these two phenomena give 
rise to unprecedented non-reciprocal responses in the context 
of LWAs. We note that relatively similar approaches have also 
been recently applied to develop non-reciprocal LWAs at 
microwaves [62], [173]. 


Such exciting responses can be achieved at THz simply 
using an antenna identical to the one of Fig. 9, but now making 
the bias signal in each pad oscillate at a frequency fın, which 
creates different conductivity profiles along the direction y at 
different instants, i.e. a spatiotemporal modulation of 
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Fig. 12. Non-reciprocal graphene-based leaky-wave antenna based on 
a spatiotemporally modulated surface impedance [43], [46]. (a) 3D 
schematic of the structure. The antenna is like the one described in Fig. 
9a, but now each gating pad is biased with a time-varying signal. The 
inset on the right illustrates the resulting graphene’s conductivity profile, 
with a spatial period p and a modulation frequency fm, at instants to, 
ti=to+1/(4fm) and tz=to+1/(2fm). (b) Dispersion relation of the n=-1 
harmonic supported by the antenna in transmission [propagation towards 
+y in panel (a)] and reception [towards -y]. Black dashed lines show the 
light cone. (c) Antenna gain in transmission and reception. Insets show 
the frequency conversion involved in both operation modes. Parameters 
are fo=350 GHz, fw=40 GHz, graphene relaxation time T = 0.5 ps and 
unmodulated chemical potential 4.9 = 0.3 eV. The structure is in free- 
space and the radiation angle @ is measured from broadside. © 2016 
IEEE. Reprinted with permission, from [43], [46]. 


graphene conductivity, as shown in the inset of Fig. 12a [43]. 
This creates and effective moving medium for the propagating 
plasmons able to break time reversal symmetry (reciprocity), 
enabling for instance transmit-only or receive-only antennas. 
Keeping the same notation as in the static case, the surface 
reactance can now be expressed as 


X's = Xz [1 + Msin (= y- 2nfast)]. (7) 


Importantly, this time-modulated surface reactance can easily 
be achieved at THz in moderately doped graphene (|,| > 
kgT), a scenario where the modulation of the chemical 
potential uo linearly translates to a similar modulation of 
graphene’s conductivity [43], [174]. In addition, graphene’s 
conductivity can follow the gating signal variations up to 
hundreds of gigahertz [73], guaranteeing an adequate 
implementation of this approach. The dispersion relation of 
the modulated plasmons can be found analytically [43], [162], 
[164]. Fig. 12b shows the radiating harmonics for a 
spatiotemporally modulated case (blue, f,, > 0) and for static 
modulation (dashed red, fm = 0). When fm > 0, the radiating 
harmonic is upshifted (downshifted) in frequency by an 
amount equal to fm for propagation towards +y (—y), 
implying different radiation properties in transmission and 
reception. Indeed, the result of the spatiotemporal modulation 
is an asymmetry in the dispersion diagram with respect to the 
frequency axis, which is a clear demonstration of 


nonreciprocal response. This exciting response can be better 
understood through Fig. 12c, which shows the gain of this 
LWA upon time-reversal, demonstrating up to 20 dB of 
contrast in transmission/reception for realistic graphene 
parameters, isolating the antenna. We stress again that the 
nonreciprocal response achieved is twofold: first, the gain and 
radiation pattern are significantly different in transmission and 
reception modes, and second, the SPPs transmitted/received 
oscillates at different frequencies under time reversal. Note 
that this non-reciprocal LWA maintains the beam scanning 
capabilities of its reciprocal counterpart, and it can 
electronically switch between standard reciprocal operation, 
transmit-only, and receive-only modes. 


Given the recent interest by the engineering community in 
non-reciprocal devices and the ease of applying space-time 
modulations to graphene, we envision the development of 
exciting plasmonic antennas and guided devices based on this 
paradigm in the coming years. 


V. REFLECTARRAY ANTENNAS 


Reflectarrays are a very well stablished topic in the RF, 
microwave, and antenna research community [175]-[177], 
and have enabled many terrestrial and satellite communication 
systems with outstanding capacity gain in the last decades 
[178]-[180]. These structures comprise an illuminating feed 
antenna and an array of reflective unit cells that each introduce 
a certain phase-shift upon reflection of a wave on the surface. 
The general concept is depicted in Fig. 13a [175]. The 
combined effect of all cells allows to engineer the phase fronts 
of the reflected beam, thus synthesizing highly directive far 
field beams in nearly any direction. It is worth noting that both 
reflectarrays and transmitarrays [181] (where the phase-shift 
is applied in transmission rather than in reflection [48]) are 
based on the same phenomenon as the so-called “generalized 
Snell laws’ that have inspired abundant research by the optics 
community in the past few years [181]-[186]. Reflectarrays 
combine the advantages of parabolic reflectors and phased 
arrays, forming narrow, reconfigurable beams through planar 
light-weight structures with simpler feeding networks than 
usual array antennas. Different technologies have been 
proposed to dynamically control the phase of the unit cell at 
micro and millimeter waves, such as semiconductor diodes 
[187]-[189], MEMS lumped elements [180], [190], [191], or 
liquid crystals [192], [193]. They are, however, not suitable 
for the THz band mainly due to size and loss. Graphene is a 
very promising candidate to fill this void thanks to its high 
inductance and strong electric field effect that provide 
miniaturization and the required tunability to synthesize any 
reflection phase. 


The first graphene-based reflectarray, illustrated in Fig. 
13b, was introduced by E. Carrasco et al [47], [48]. The 
proposed unit-cell is composed of a square graphene patch 
transferred onto a S1O2 substrate, as further described in Fig. 
13b. In conventional reflectarray patches made of gold or 
noble metals, their resonance occurs when the size of the patch 
is comparable to a half wavelength of the wave in the effective 
media, which usually leads to large structures. This scenario 
is different in the case of graphene patches, since, as 
previously discussed in Section I, their resonant frequency is 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


+f. 


Graphene 
y Z Alumina 


= "> Polycrystalline Silicon 


Source 


— Quartz 


= = p= 006V } 
sensns fle = 0.29 eV 






Gain (dBi) 








-30-20-10 0 10 20 30 40 50 60 
Elevation (deg) 


(d) 


Fig. 13. Reconfigurable graphene-based reflectarray [48]. (a) 3D 
schematic of the concept. (b) Details of a unit-cell element. (c) Phase 
shift (in degrees) induced in reflection at 1.3 THz by square patches of 
different sizes versus graphene chemical potential. (d) Antenna radiation 
diagram (E plane) at 1.3 THz upon reconfiguration. Graphene relaxation 
time is 1 ps and other parameters are a=b=14 jum, a,=b, in the range of 
0-10 jim, and quartz, alumina, and polysilicon thickness are 30, 0.1 and 
0.05 um, respectively. © Reprinted from [48], with the permission of 
AIP publishing. 


much lower than their metallic counterpart due to the 
plasmonic response of graphene. For instance, in the example 
of [48], the size of the required graphene patches is below 
)/10, allowing to achieve an impressive miniaturization factor 
of around 5. Furthermore, a DC voltage applied between the 
graphene patch and the polycrystalline silicon layer allows to 
tune the patch reactance and thus the reflection phase of the 
unit cell. We stress that a thorough knowledge of the unit cell 
phase response, both versus its size or tunable features, is the 
most important step in any reflectarray design [175]. Fig. 13c 
shows the reflection phase (in degrees) for the tunable 
graphene patch versus chemical potential and patch size. 
Interestingly, the dynamic range of phase variation is much 
larger around low values of ue due to the dependence of the 
plasmon wavenumber on uç, making them better suited for 
reflectarray design. This is qualitatively explained by Eq. 1, 
which shows that for large uw, the plasmon wavenumber 
asymptotically approaches ky. Once the reflection properties 
of the unit cell versus the design variables have been 
determined, the phase distribution that implements the desired 
far-field pattern can be readily found [47], [48]. In this case, 
the final reflectarray can be designed by varying both 
parameters (patch size and wu.) or keeping either of them 
constant. Fig. 13d shows the radiation pattern of a reflectarray 
based on size-modulation for u, = 0.0 eV and p, = 0.29 eV 
illustrating the excellent directivity and beam-scanning 
capabilities of this technology in the terahertz band. 


One exciting application of reflectarrays is the generation 
of vortex beams [194]-[200]. Vortex beams are characterized 
by their orbital angular momentum (OAM), and they allow to 
increase the spectral efficiency thanks to the orthogonality of 
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Fig. 14. Reconfigurable graphene-based reflectarray for vortex beam 
generation [206]. (a) 3D schematic of the structure. Dark and bright gray 
colors correspond to unit cells with different chemical potentials. (b) 
Simulated 3D radiation patterns (left) and phase front (right) of 
generated vortex beams with modes /=0 (top) and /=-1 (bottom). Details 
of the geometry and applied voltages can be found in [206]. © 2016 
IEEE. Reprinted with permission, from [206]. 





their ideally infinite number of states. Vortex beams have 
attracted significant attention in the past years thanks to 
promising experiments claiming that OAM might provide 
enhanced transmission capacity than any other known 
communication system [201]. Recently, it was clarified that 
OAM multiplexing is in fact a subset of MIMO systems [202], 
[203], and consequently, cannot bring larger capacity gain that 
them. Nevertheless, vortex beams can lead to exciting 
applications [204], for instance in communications with 
enhanced encryption capabilities, optical tweezers, optical 
trapping, LADAR (laser detection and ranging), and they 
propagate with less distortion than common Gaussian beams 
through turbulences. Unfortunately, the generation of vortex 
beams is a challenging task [196], [197], [205]. A common 
approach based on reflectarrays technology takes advantage 
of the reflection from a plane wave that impinges over a 
sectored circular surface, with each sector having the same 
reflection magnitude while their phase increments add up to a 
multiple / of 27 (/ determines the order of the vortex beam). 
The required reflection coefficient of the nth sector and the 
corresponding surface impedance can be found analytically 
[206]. 

Since vortex beam generation requires precise control of the 
reflection phase, graphene patches offer once again a perfectly 
suited platform with large design flexibility. This concept is 
illustrated in Fig. 14 [199]. The structure is qualitatively like 
the previous reflectarray, but the patch parameters now divide 
the circular surface in phase sectors. Like before, the key step 
in the design is the characterization of reflection properties 
versus the unit cell parameters, which follows the same 
guidelines as discussed above and is omitted here for brevity 
[206]. In the pictorial representation of Fig. 14a, each shade 
of gray corresponds to patches with different u,. Fig. 14b 
shows the radiation pattern and phase fronts of the generated 
vortex beams with modes | — 0 and l = —1, but we note that 
higher order modes can also be implemented by simply 
changing the biasing scheme. Importantly, this structure can 
also function as a regular beam-scanning reflectarray using 
appropriate bias values. 


The same reflectarray concepts can be used to steer 
reflected or transmitted beams at will without necessarily 
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Fig. 15. Reconfigurable reflectarray based on aperiodic array of 
graphene ribbons [49]. (a) 3D schematic of the structure. (b) Phase of 
the reflection coefficient as a function of graphene chemical potential 
for different widths W of the nanoribbons. Mutual coupling between 
identical neighboring ribbons is considered. (c) Width of the 
nanoribbons across the device. (d) Far-field radiation pattern for 4 = 
1.0eV_ when a Gaussian beam illuminates the device at an incident 
angle of 45. The reflected beam can be steered by tuning 4. Further 
details of the geometry can be found in [49]. © 2015 IOP Publishing. 
Reprinted with permission, [49]. 
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implementing an antenna, akin to optics papers that exploit the 
generalized Snell’s laws [182]. Fig. 15 implements light 
bending in reflection by means of an aperiodic array of 
graphene nanoribbons [49]. To design such device, it is 
essential to fully understand the reflection phase of individual 
strips, both as a function of its width and chemical potential, 
as investigated in Fig. 15b. Using this information, an 
appropriate selection of the strip widths permits to implement 
a linear phase discontinuity profile able to reflect an incident 
beam to the desired direction. Fig. 15c-d show the strip width 
distribution of the final reflectarray and a far-field radiation 
diagram, respectively. This configuration also provides 
tunable features, allowing to direct the reflected beam to 
almost any direction in the space by using a unique DC bias, 
thus greatly simplifying the complex biasing scheme of 
previous graphene-based reflectarray designs [47], [48]. 
Lastly, we stress that all these concepts can also be 
implemented in transmission through transmitarrays [181], 
[207], but they are usually less suited for antenna applications 
due to increased difficulty in fabricating the bias circuitry and 
the requirement of using multilayers to achieve wide phase 
control with minimal reflections [181]. 


VI. FEEDING NETWORKS 


Previous sections have exclusively discussed graphene- 
based antennas, i.e. structures able to couple guided energy to 
far field radiation (or vice versa). However, future THz 
transceiver front ends will also require guiding components 
such as interconnects, modulators, switches, filters, and phase 
shifters with high degrees of tunability enabled by graphene 
to properly feed the radiating elements in a unique, 
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Fig. 16. Graphene-based components for THz antenna feeding 
networks. (a) Plasmonic switch based on graphene's field effect [213]. 
(b) Low-pass THz filter [97]. (c) Band-pass filter [214]. (d) Phase- 
shifters [215], [214]. Inset on the top of panels (b) and (d) describe 
potential implementations of the devices. Further details of the 
geometries can be found in the literature. Panel (a): O 2013 OSA. 
Reprinted with permission [213]. Panel (b): O 2014 IEEE. Reprinted 
with permission, from [97]. Panels (c), (d): 2014 American Chemical 
Society. Reprinted with permission, from [214]. 


miniaturized, and multifunctional device. Clearly, an adequate 
interconnection between graphene devices and antennas will 
lead to lower losses and mismatch that using alternative 
technologies, plus it might also provide enhanced responses 
thanks to graphene's exotic properties. In this section, we 
review some of the most relevant advances in this area. The 
components described in the following exploit in some way 
the influence of graphene's chemical potential on the 
propagation constant and characteristic impedance of SPPs in 
graphene waveguides (see Eqs. 1 and 2), which allows to 
implement a very wide range of THz transmission lines by 
simply changing a biasing signal. This simple concept has 
tremendous implications, as virtually any textbook device that 
can be modeled using a transmission line approach can 
therefore be implemented using graphene plasmons in a 
miniaturized and reconfigurable THz platform [1], [25], 
[208]-[212]. 

Fig. 16 shows promising examples of several fundamental 
guided components, namely a switch [213], low-pass filter 
[94], band-pass filter [214] and a phase shifter [215], that are 
all based on the propagation of SPPs along graphene 
waveguides. The switch, illustrated in Fig. 16a, operates by 
introducing a strong impedance mismatch in a short graphene 
section with different chemical potential, easily obtained by 
biasing a gating pad located beneath that region. This simple 
approach allows to simultaneously achieve large isolation 
(240 dB), low insertion loss, and fast switching speed thanks 
to graphene's field effect [73], [213], [216]. The low-pass and 
bandpass filters, shown in Fig. 16b-c, use a similar principle 
to implement traditional microwave filtering functions, where 
specific characteristic impedances and electrical lengths must 
be implemented at the desired frequency [214], [217], [218]. 
Specifically, each of the gated graphene sections can be 
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Fig. 17. Influence of spatial dispersion in the response of graphene- 
based phase shifters (top row) and low-pass filters (bottom row) [224]. 
The analysis of the phase shifter includes (a) phase difference between 
ports and (b) phase error due to spatial dispersion versus the permittivity 
of the surrounding media. The 7th-degree low pass filter is designed and 
analyzed in (c) free space and (d) embedded in Si (€, ~ 11.8). Solid 
lines: results neglecting spatial dispersion effects; dashed lines: results 
including spatial dispersion effects. Further details of the devices can be 
found in [224]. O 2014 IEEE. Reprinted with permission, from [224]. 


modeled as a tunable transmission line, thus allowing to fully 
exploit decades of research in the 20" century dedicated to 
refining mathematical functions for optimal microwave filters 
[218]-[220]. Remarkably low insertion loss and band 
rejection are achieved in both structures, which also exhibit 
unprecedented electronic tuning of the filter central frequency. 
Fig. 16d illustrates a plasmonic phase shifter able to 
implement a wide range of phase shifts while providing small 
return loss, insertion loss, and phase error. These are just 
pioneer examples of potential guided devices, and that we 
expect many more microwave and optics devices to be 
successfully translated into THz graphene plasmonics soon. 


Usual analysis and design of graphene-based THz antennas 
and guided components have employed standard techniques 
borrowed from the microwave and optics communities [1], 
[2], [221]-[223]. However, there is a key difference when 
dealing with components based on ultra-confined surface 
plasmons polaritons (i.e., Kspp >> ko, where ks;pp and Kg are 
the plasmon and free-space wavenumber, respectively) 
propagating in graphene: the influence of nonlocal responses 
or spatial dispersion [82], [96], [116], [209], [224]-[230]. 
This phenomenon is related to the finite Fermi velocity of 
electrons moving through the material and it fundamentally 
limits the maximum plasmon confinement. It arises because 
carriers on graphene have a fixed maximum speed (Vp =% 
c/300), and therefore they cannot follow sharp spatial 
variation of ultra-confined surface plasmons, limiting the 
maximum SPPs wavenumber to be kspp < 300k, at any 
given frequency. In fact, spatial dispersion is non-negligible 
even for lower values, being significant in graphene when 
Kkspp = 100K . We do stress that the spatial dispersion 
mentioned here — which effectively translates onto a 
wavenumber-dependent graphene conductivity 
a(Kspp) [109], [110] — is intrinsic to the band structure and 
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Fig. 18. Magnet-less non-reciprocal graphene guided device [43], 
[174]. (a) Schematic of a spatiotemporally modulated graphene 
waveguide. The conductivity of the bottom graphene layer is modulated 
through the time-varying voltages V;(t). The inset illustrates the 
resulting graphene's conductivity profile at different temporal instants. 
(b) Dispersion diagram of the unperturbed fundamental modes in a 
graphene-based parallel plate waveguide. Arrows indicate the frequency 
and wavenumber shifts imparted by the spatiotemporal modulation. (c)- 
(d) Electric field (longitudinal component) in the isolator at frequencies 
fo. fo + fm when the device is excited at f; from (b) the left and (c) the 
right sides (M is the modulation depth). Further details of the device can 
be found in [43]. € 2016 IEEE. Reprinted with permission, from [43], 
[174]. 









properties of graphene, and should not be confused with the 
spatial dispersion that might arise in metamaterials due to the 
finite size of their composing unit-cells. Importantly, 
nonlocality might affect some devices more than others, 
depending on their features and operation principle. In the 
case of THz antennas, their radiating capabilities are always 
linked to their electric size in terms of free-space wavelength, 
which forces trade-offs between miniaturization and antenna 
figures of merit like gain. For this reason, ultra-confined SPPs 
are usually avoided in antennas, since their coupling to far- 
free radiation can be challenging. This is in principle not the 
case for guided devices. Their dimensions can eventually be 
reduced as much as desired, as long as in- and out-coupling of 
power can be done efficiently, so ultra-confined SPPs find a 
key role in these devices. Such SPPs can easily be obtained by 
decreasing the chemical potential, increasing the substrate 
permittivity (see Eq. (1)) of graphene waveguides, or using 
extremely anisotropic and hyperbolic configurations [231]- 
[235]. Thus, incorrectly assuming a local graphene response 
can lead to very significant errors in the design of guided 
component such as phase shifters and filters, as illustrated in 
Fig. 17 [224]. For instance, Fig. 17a-b show the phase error 
that may appear in graphene phase shifters due to spatial 
dispersion, demonstrating a strong dependence with the 
substrate used. Fig. 17c-d consider the influence of spatial 
dispersion in low-pass filters like those of Fig.16. In the free- 
standing case (Fig. 17c), the local conductivity model is a very 
good approximation, whereas performance is severely 
affected in Fig. 17d, which uses a high permittivity substrate. 
The operation frequency is shifted, return loss increased, roll- 
off decreased, and the implementation requires higher 
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voltages. Similar issues are expected for extremely 
subwavelength plasmonic switches [213], [236]-[238] and 
other guided devices. In general, spatial dispersion limits the 
maximum tunable range of graphene plasmonic devices, 
reduces their quality factor, and degrades their frequency 
response [96], [116], [209], [224]. However, these effects can 
usually be considered and counteracted in the design of 
graphene-based devices by using accurate nonlocal models of 
graphene conductivity [116], [224]. 


Non-reciprocal waveguides are also of great importance for 
future THz devices to isolate sources from unwanted 
reflections or to enable full-duplex communications. Several 
groups have proposed non-reciprocal components based on 
magnetically biased graphene [52], [54], [111], [239]-[241], 
which may from significant drawbacks, namely the potential 
large size of the required magnets and their lack of 
compatibility with integrated circuits — in practical 
applications. A more promising solution for feeding networks 
is to adopt the magnet-less spatiotemporal modulation scheme 
described in Section III and applied there to develop non- 
reciprocal LWAs [43], [46]. This strategy might also allow to 
achieve non-reciprocal responses in both radiating elements 
and feeding networks, thus increasing isolation to almost any 
desired level. Fig. 18 illustrates the operation of a 
spatiotemporally modulated graphene parallel plate 
waveguide that implements a waveguide plasmonic isolator 
[43], [82]. By appropriately modulating the conductivity of 
one of the graphene layers in space and time, unidirectional 
mode conversion can be realized. Fig. 18b depicts the 
operation principle. It shows the dispersion of the two 
fundamental modes of the unmodulated structure, as well as 
the non-reciprocal wavenumber and frequency shifts induced 
by the modulation (black arrows). The modulation enables 
phase-matching and thus full exchange of energy in the 
backward direction (propagation towards -z), but not in the 
forward one (propagation towards +z). This isolating 
phenomenon is visually apparent in Fig. 18c-d. When the 
structure is excited from the left (towards +z), the modulation 
does not phase-match the modes, so only very small residual 
coupling occurs. When excited from the right (towards -z), 
phase-matching is perfect, and all power at the original 
frequency fo is converted to a different mode and frequency, 
isolating the ports. The residual power at fp + fj, is entirely 
coupled to an orthogonal mode that can be readily scattered or 
reflected with modal filters. We stress that this pioneering 
example is just the first magnet-less non-reciprocal guided- 
device using THz graphene plasmonics, and we expect more 
refined configurations to be developed, such as integrated 
circulators, as well as theoretical studies on the fundamental 
limits of non-reciprocal responses. 


VII. APPLICATIONS 


The most immediate application for graphene antennas is in 
THz transmitter front ends for communication systems. 
Successfully combining graphene antennas with passive 
plasmonic devices and efficient THz sources will lead to 
unprecedented wave manipulation functionalities in the THz 
band, enabling picocell and interchip wired or wireless 
communications with ultra-high bandwidth [16], [18], [19]. 
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Fig. 19. Conceptual nano-transceiver architecture (a) and a specific 
implementation of a plasmonic nano-transmitter (b) based on high 
electron mobility transistor (HEMT) built with a III-V semiconductor 
and enhanced with graphene [242]. O 2014 IEEE. Reprinted with 
permission, from [242]. 


Although efficient excitation and detection of THz plasmons 
is still a challenging problem, significant strides have been 
made in parallel to the development of the passive devices that 
occupy this review. Excitation schemes used by most groups 
around the world do not currently employ graphene, using 
instead variations of conventional near-field techniques that 
have been shown to excite plasmons up to a hundred times 
shorter than the free-space wavelength [63], [64]. 
Nonetheless, some groups have proposed graphene-based 
sources that aim to maximize compatibility with graphene 
feeding networks and antennas with minimal interconnection 
[140], [242]. Fig. 19 shows an example of this approach based 
on a High Electron Mobility Transistor (HEMT) [242]. In 
transmission, the applied voltage between the drain and source 
accelerates a beam of electrons in the HEMT channel that 
excite a plasma wave able to couple to graphene SPPs. In 
reception, reciprocity allows to simply measure a drain-source 
voltage, more easily detected than near fields, especially in 
integrated systems. The graphene SPPs generated in this way 
can readily drive resonant and leaky-wave antennas, or any 
other passive device with good impedance matching. 


Besides THz communications, one of the most promising 
applications of graphene plasmons is in biological sensing, 
since DNA molecules and many large proteins have rotational 
and vibrational modes in the THz and far-IR [17], [214], 
[243]-[249]. Plasmonic sensing has a long history at optical 
frequencies using noble metals like gold and silver [160], 
[250]-[254], but this technology could not be translated to 
THz before the advent of graphene. Moreover, the thinness 
and optical transparency of graphene make it very attractive 
for wearable and implantable sensors. The extreme field 
enhancement enabled by graphene plasmons in the THz band 
makes electromagnetic observables such as transmittance, 
phase, or polarization extremely sensitive to small 
perturbations in the complex permittivity of the surroundings 
[255], [256]. This sensitivity can be exploited to detect small 
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Fig. 20. Graphene-based THz bio-sensors. (a) Schematics of a THz 
biosensor composed of a nanofluidic delivery channel and the structure 
shown in Fig. 16c [214]. (b) Schematics of a graphene-based harmonic 
sensor designed for various point-of-care monitoring and wireless 
biomedical sensing [257]. The right panel illustrates an eye-wearable 
device (smart contact lens) based on the device, which may detect in real 
time the pathogen, bacteria, glucose, and infectious keratitis. Panel (a): 
© 2014 American Chemical Society. Reprinted with permission, from 
[214]. Panel (b) Reprinted from [257], with permission of AIP 
Publishing. 


concentrations of particles that respond strongly to THz light. 
Fig. 20a depicts a THz biosensor concept based in this 
principle, where the biological agents to be analyzed flow 
through a graphene plasmonic waveguide [214]. The presence 
of these molecules causes small changes in the channel’s 
properties, detectable through the scattering parameters of the 
waveguide. This is especially suitable for molecules with THz 
absorption resonances, as they have a strong effect on the 
channel’s permittivity. Note that, in this sensor, SPPs in the 
waveguide must be launched by either an external near-field 
source or an antenna. In contrast, Fig. 20b shows a recently 
proposed device where all the required elements (the sensor, 
frequency modulator, and antenna for energy harvesting) are 
included in a single module [257]. The structure is designed 
as a dual resonant antenna part of a graphene field effect 
transistor (GFET). The rectifying mechanism of the GFET 
generates a second harmonic signal whose intensity depends 
on the concentration of molecules to be sensed (this is a very 
sensitive process [257]). Since the monopole antenna is 
designed to resonate at both the original and the second 
harmonic frequency, the sensing output is automatically 
reradiated and available for processing. 


We envision future research on THz sensors to follow 
several directions: first, using metasurface concepts like 
hyperbolic dispersion to further enhance the near fields for 
increased sensitivity [117], [231], [232], [235], [258]-[264], 
second, designing and optimizing graphene nanoparticles 
supporting  ultra-sensitive localized surface plasmon 
resonances (the most sensitive plasmonic nano-sensors at 
optics are based on nanoparticles [250]), third, using nonlinear 
responses [265]-[273], and fourth, employing nano- 
mechanical resonators [274], [275]. All of this must be 
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Fig. 21. Graphene-based THz detector based on a log-periodic 
circular-toothed antenna [70]. (a) 3D schematic of the structure. Off-axis 
parabolic mirrors focus the THz beam into the device. The detector 
consists of a log-periodic circular-toothed antenna patterned between the 
source and gate of a GFET. The drain is a metal line running to the 
bonding pad (b) Photoresponse signal of the detectors a function of 
incident terahertz power. Dots are the experimental data; the dashed line 
is a fit to the data. Adapted by permission from Macmillan Publishers 
Ltd: Nature Materials [70], copyright 2012. 


accompanied by suitable instrumentation capable of high- 
throughput detection. In this context, graphene THz detectors 
have been subject of study for a few years, and promising 
advances have been made theoretically and experimentally. 
Current commercially available THz detectors are either 
based on thermal sensing elements or on nonlinear electronics. 
The former are typically slow or require cryogenic cooling, 
while the latter are limited to frequencies lower than 1 THz 
[276]-[278] . A more recent approach relies instead on field- 
effect transistors where a plasma wave is excited in the 
transistor channel by the impinging electric field, potentially 
enabling faster and more sensitive detection [279]-[283]. 
Different technologies have been used as platform for this 
class of FETs, including II-V  high-electron-mobility 
transistors, Si-based MOS, and InAs nanowires [276], [278]. 
However, their performance is still inferior to other 
technologies, and is limited to the low THz. Graphene holds 
the promise to overcome these issues due to its higher electron 
mobility and suitability for manipulating THz waves. A recent 
example working at room temperature is depicted in Fig. 21 
[70]. A log-periodic circular-toothed is patterned between the 
source and gate of a GFET. The concept is similar to the 
single-module sensor discussed above and the THz 
transceiver from Fig. 19, but the structure is now optimized 
for operation as a detector. When a THz oscillating electric 
field is fed between the gate and the channel of the GFET, the 
wave is rectified by the nonlinearity of the transfer function 
[70], inducing a DC signal between source and drain that is 
proportional to the received optical power. Fig. 21b plots the 
photoresponse versus incident THz power in fabricated 
GFETS, demonstrating voltages of the order of microvolts for 
incident powers of 1 miliwatt. The NEP of this detector is 
reported around «30 nW Hz "7, and the authors also 
demonstrate THz imaging with daily life objects [70] — THz 
waves are ideal for imaging common dielectric substances, 
hence their importance in security-related applications [17]. 
Terahertz imaging also has promise for label-free DNA 
analysis methods with increased diagnosis accuracy and lower 
cost, and in medical diagnosis for detection of anomalous 
tissue [17], [284]-[286]. There is still much room for 
improvement in GFET detectors, most notably by perfecting 
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Fig. 22. Graphene-based THz modulators [72], [76]. (a) Conical band 
structure of graphene and associated optical processes. Intraband 
transitions (red arrows) dominate under terahertz illumination [72]. (b) 
Schematic of a single-layer graphene terahertz modulator. A monolayer 
graphene sheet was transferred on a SiO2/p-Si substrate [72]. (c) 
Measured (symbols) and modelled (lines) intensity transmittance for a 
single-layer modulator as a function of frequency for back gate voltages 
of 0 V (blue) and 50 V (red) [72]. (d) Schematic of a double-layer 
graphene modulator, where the layers gate each other [76]. Inset (top 
right) shows an illustration of the fabricated device. (e) Measured 
conductivity of the stack plotted versus a voltage Vpc applied between 
the two graphene sheets [76]. Simulated results are included for 
comparison. Panels (a)-(c) adapted by permission from Macmillan 
Publishers Ltd: Nature Communications [72], copyright 2012. Panels 
(d),(e) adapted by permission from Macmillan Publishers Ltd: Nature 
Communications [76], copyright 2015 


graphene fabrication processes, but the outlook is certainly 
promising, and this technology could soon have a real impact 
on society [71], [287]-[291]. For instance, sandwiching 
graphene between hexagonal boron nitride layers leads to 
significant improvements in mobility as has been recently 
demonstrated by several groups in other contexts [63], [68]. 


The last major class of device we will consider, the 
modulator, will also be essential for a transition towards 
practical THz high speed communications. THz modulator 
technology is in its infancy, and until recently their 
development was lagging behind passive devices and 
sources/detectors [292]. Many attempts have been made at 
semiconductor modulators, but they always suffered from 
either low modulation frequency or low modulation depth, 
limited by the achievable tunability in carrier density. 
Graphene is the ideal material to overcome these two issues, 
since carriers can be injected very fast, holes conduct as well 
as electrons, and it interacts strongly with THz fields [26], 
[73]. As one may expect, the operation of graphene 
modulators is usually based on rapidly tuning the chemical 
potential using a dynamic electronic bias to modulate the 
conductivity and thus the transmission. Pioneering 
experiments by Sensale-Rodriguez et al demonstrated a few 
years ago substantial improvements over the state of the art in 
a simple modulator [72], [77], [293]. The structure, depicted 
in Fig. 22a-c, comprises a graphene monolayer transferred on 
a SiO;/p-Si substrate and biased with metal gates. A voltage 


applied between the gates injects carriers on graphene, 
modulating the chemical potential. Fig. 22c plots the 
measured intensity transmittance as a function of frequency 
for back-gate voltages of O0 and 50 V, showing a 1596 
modulation depth at around 570 GHz. This performance was 
remarkable at the time, especially considering the simplicity 
of the design and the poor-quality graphene used in the 
experiment. A promising alternative to metal gates relies 
instead on self-biased graphene stacks where the tuning 
voltages are applied directly between graphene layers, as 
illustrated in Fig. 22d (inset shows the fabricated sample) [76], 
[127]. This approach minimizes the footprint and potentially 
reduces the RC constant of the structure, enabling faster and 
smaller devices. For small separation between the graphene 
sheets, the structure behaves as a single ultra-thin layer with 
conductivity equal sum of the two - this concept can be 
extended to an arbitrary number of layer pairs [72], [76]. The 
reconfiguration capabilities of the stack are plotted in Fig. 22e, 
demonstrating ambipolar tunability of the overall stack 
conductivity. Incorporating metal gates or additional graphene 
sheets to the modulator further boosts the available range of 
conductivity configurations, as discussed in [76]. Other 
modulators have been proposed and fabricated at THz using 
the same principles [74], [99], [292], and also at infrared 
frequencies using integrated waveguides [73], [75]. At 
infrared, graphene offers similar tunability and it can be easily 
integrated with the low-loss, technologically mature silicon 
photonics structures. In Ref. [73], the authors demonstrated 
modulation speeds of up to 30 GHz in hybrid dielectric- 
graphene waveguide modulators [73] using low quality 
graphene — improvements in graphene quality and low- 
resistance contacts may allow speeds in the hundreds of GHz 
[73]. Importantly, graphene's loss is less critical in this type 
of structure because its function is to perturb the photonic 
modes (with very high quality factor), rather than to be the 
main mechanism of energy propagation. These results are also 
a very strong indicator of the far-reaching, untapped potential 
of this hybrid approach, which we predict will be exploited in 
coming years not only to implement modulators, but also as 
an enabler of device tunability and non-reciprocal responses 
via spatiotemporal modulation in low-loss photonic 
technology. 


VIII. CONCLUSION AND OUTLOOK 


We have reviewed the evolution of graphene antennas, from 
the inception of this field to the present state of the art, 
discussing their operation, shortcomings, and prospects. Two 
common themes stand out throughout this review: 
reconfigurability and miniaturization. These are the most 
unique and attractive features of graphene, and they are the 
cornerstone that has allowed researchers to conceive widely 
tunable antennas, ultra-fast modulators and detectors, 
sensitive plasmonic sensors, and magnet-less non-reciprocal 
plasmons. Although the practical realization of most of these 
devices has so far proved difficult, recent progress in graphene 
fabrication processes suggests a promising future. This is 
without a doubt the most important missing piece towards 
feasible graphene technology. From the conceptual point of 
view, we believe research will follow several directions in the 
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coming years. First, we expect an increased cross-fertilization 
with the field of metamaterials and metasurfaces, exploiting 
the exotic electromagnetic phenomena that they enable. A 
second likely direction consists of using graphene’s ultra-fast 
reconfigurability to introduce spatiotemporal modulations 
able to break time-reversal symmetry, enabling a wide variety 
of non-reciprocal devices. A third promising direction consists 
of combining graphene plasmonics with other technologies 
like silicon photonics to reduce the sometimes-stringent 
requirements on graphene quality. We envision, in the long 
term, that all these pieces will converge into a fully integrated, 
miniaturized, multifunctional, and reconfigurable graphene 
platform able not only to push forward but to transform 
current boundaries of THz technology. 
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